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УДК 539.3 


Трёхмерная контактная задача о взаимодействии упругого клина 
с двумя штампами при учёте трения 


Д. А. Пожарский, А. С. Поляков 
(Донской государственный технический университет) 


В квазистатической постановке изучена пространственная контактная задача для упругого клина, в одну 
грань которого симметрично вдавливаются два одинаковых эллиптических в плане жёстких штампа с учётом 
трения. Другая грань клина свободна от напряжений. Область контакта считается неизвестной. Рассмотрены 
случаи движения штампов перпендикулярно ребру клина. При этом штампы могут приближаться к ребру или 
удаляться от него. При использовании известной функции Грина для трёхмерного клина контактная задача 
сведена к интегральному уравнению (ИУ) относительно контактного давления. Для решения ИУ применяется 
метод Галанова, позволяющий одновременно определить область контакта и давление в этой области. Сде- 
ланы расчёты при разных коэффициентах трения, углах раствора клина и расстояниях между штампами. 
Ключевые слова: теория упругости, контактная задача, клин, трение, метод Галанова. 


Введение. Исследована задача о квазистатическом контакте с трением двух симметрично распо- 
ложенных штампов на грани трёхмерного упругого клина. Другая грань клина свободна от 
напряжений. Штампы начинают медленно двигаться перпендикулярно ребру клина. Симметрич- 
ные области контакта неизвестны. Для решения использован метод нелинейных граничных ИУ [1, 
2]. Ранее исследовались аналогичные задачи для одного штампа без учёта трения [1] и с учётом 
трения [3]. Изучалось взаимодействие двух штампов без учёта трения [4, 5]. 

Постановка задачи. Рассмотрим в квазистатической постановке трёхмерную контактную задачу 
для упругого клина {ле [0, ©]; фе [-а, а]; 2 Е (-00, со)} с учётом сил трения, действующих пер- 
пендикулярно ребру клина, когда грань клина ф = —а свободна от напряжений. Пусть два сим- 
метричных жёстких штампа, первоначально внедрённых в грань ф = а, начинают достаточно мед- 
ленно двигаться по этой грани (без перекоса) перпендикулярно ребру клина. Упругие параметры 
материала клина — С (модуль сдвига) и у (коэффициент Пуассона). Силы трения коллинеарны 
направлению движения (и направлены против движения). Штампы имеют форму эллиптических 
параболоидов, их подошвы описываются функциями 


(г -а) (2+ в} 

2, ы 2К 
Задача симметрична относительно полуоси г. На штампы действуют нормальные вдавливающие 
силы Р, имеющие плечи [Н. (относительно полуоси г) и Н, (относительно ребра клина), и касатель- 
ные силы 7, направленные перпендикулярно ребру клина. В соответствии с законом Кулона вы- 
полняется соотношение 7 = |рР, где и — коэффициент трения Кулона. Условие контакта штампов 
и клина имеет вид: 





Г [уе = р В (1) 


2 


и, (га, 2) =-[6-Е, (г,2)], (г,2)е®,, (2) 


где ц,(г, а, 2) — нормальное перемещение упругого материала клина с одной свободной от 
напряжений гранью при действии на другой грани нормальных и касательных напряжений, м; 
б — осадка штампов, м; ©. — неизвестные области контакта, м. Смена знака и эквивалентна 
изменению направления движения штампа на противоположное. Трением вдоль ребра клина 
пренебрегаем. 
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При известных величинах а, |, С, \, д и заданных функциях Ё(г, 2) требуется определить 
нормальное контактное давление оз(/, а, 2) = —-4(/ 2), (/, 2) е ©., найти области контакта ©., а 
также значения Р, Н;и Н.. 

Решение задачи. При помощи условия контакта (2) выводится ИУ относительно функции ((/, 2) 

в неизвестной области ©.. Решив это уравнение и определив ((, 2) и ©., можем найти величины 

Р, Н.и Н, из условий равновесия штампа. Для вывода ИУ рассмотрим две вспомогательные крае- 
вые задачи. Граничные условия первой имеют вид: 

ф=а: в, =-@6{г-х)5(2-у), ть=ть =0; (3) 

ф=-а: О.Т. =Т., = 0. (4) 


Здесь 0(х) — дельта-функция Дирака. Решение задачи (3), (4) известно [1]. Используя его, 
найдём 




















а - мы 
и га:2 Гаити и ‚ (Вх К. (Вг)со$В (2 - у} аиаВ, 
9=с/( ты 
5Н4аи + и зт4а си2аи = с0$2а 
и = = + 
(и) сН4аи — 2и? т? 2а-1' из (0) = наи +изт2о' (5) 
где функции Ф‚(и) (п = 1, 2) удовлетворяют ИУ Фредгольма второго рода: 
4, (т) =(1- 2) [1, в (и) + «К. (р) (О<т<®)}, (7) 
0 
Е $НпЕ д, (#) 4 
Е 2 Но но 
И "Г сиё + сппт) (сппё + сИпи) " (8) 
Са 51? 2а {ПаЁ 51? 2а 
9 (Е) > (Е) сИ2аЁ + со$4а ‘` (9) 


СН2 а — с0$4а ' 
С использованием известного интеграла [1] выделим в формуле (5) главную часть, соответству- 
ющую упругому полупространству: 


И (,а,2) =- ©, -е [ти К, (В/)| (М (и) спи) к, (Вх) + 











210 к (10) 
Ил (и) 4 (и) - И (и), (и) 
+ эсн(пи / 2) со5В (2 — у ) иаВ, 
в - (1-х) +(2-у) |. (11) 
Граничные условия второй вспомогательной краевой задачи имеют вид: 
ф=а: б,=т,, =0, т,=0@5(г-х)5(2-У); (12) 
ф=-а: Ч. =Тжь=То, =0. (13) 
Решение задачи (12), (13) также известно [3]. Используя его, найдём 
20, * и (м) ча (и) - ИА (м) (м) 
а, = = т [958 (2-у)эппик, к (Вх)Е, (и) + — ев (пи 12) диаВ м 
м К (Вг)К, (Вх) 
ее |958 (2- УВ], (ив, (Е) -и (и), (Е )]впи ЗЕ т М 
2и п? 2а (1—2\)И2аи зи п2а 
о ( И. т рая (15) 


си2аи = со$2а 
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Здесь функции 4" (и) (п = 1,2) удовлетворяют ИУ Фредгольма второго рода 


у.) @- 2) о), им (05), 








Е, (и)К, (Вх) ей, (Е)К„ (Вх) 
ое а сплести 1, 


Ва 2(1-\)(1-2\)пт2а 
(6) сп2аи = с0$2а 





Используя значение интеграла 
е СППЕЁСИПИ —1 


с [сов (2 -у) [| 


со СИПЁ- СППИ 


к, (Вг)к, (Вх) аи =-^ = Е 





выделим в (14) главную часть, соответствующую упругому полупространству: 


1-2 
о а то [еоев(а -ути К, (и; 


к, (вх) (и) и 4) ов+ 





и, (г,а,2) = 





вре |958 (2- УВ Ги (ив, (и (и), (Е) ппи эвпё - 


К, (Вг)К„ (Вх) 


-2 (1 — 2%) (сппи сНпЁ 1} сп — спи 


аи. 





(16) 


(17) 


(18) 


(19) 


(20) 


В соответствии с принципом суперпозиции заменим в (5) силу О, на искомое давление С(х, у), а 
силу © в (20) — на распределённые силы трения т(х, у) = нах, у). Проинтегрировав далее пра- 
вые части полученных соотношений по области контакта ©, ›®_ по переменным х, у и сложив 


их, получим при помощи р контакта следующее ИУ: 
Па (х,у,г,2)акду = 4, (г,2), (г,2)е9, 39, 


©. 59 


где 


[ (х,У,г,2)=Ё (х,у,г,2)+НЁ (х,У,г,2), 9, (г,2) = 210 [6 -Ё. Е 


[м (и)-сли)к, (вх) + 


Е (х,у,г,2)= 
ИИ (и)ч, (и)-и/, (и), М зори 








Е ФЫ: ый с05В (2 +у)эппиК,, (Вг)х 
0 


в (У, у, 





кои А) лв 


а |958 (2 в зав [и (ив, (Е) -и, (и), (Е) ппу 5 ПЕ — 
К, (В^)к, (Вх) 


ав. 
сппЕ — сАпи 





—2 (1 — 2%) (сппи сит -1)} 
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В силу симметрии ИУ (21) сводится к уравнению 





ау) (х,у,г,2) +, (х,у,г,2)|акау =а,(г,2), (г,2)е®.. (25) 
о 
Введём новые обозначения: 
2, =2-1, у, =У-В, а. (ге аРиЕ), оФо. (26) 
Перепишем уравнение (25) в форме 
9. (ж,у.)К (ж,уг,2,)ахау, = (г,2.), (г,2.)е ©, (27) 
о 
а} 2 
аи) Пе , (28) 
2, 2К, 
К (х,У,г,2)=Ё (х,у,г,2)+Ё, (ху +В,г, 7 +В). (29) 


Звёздочки далее опускаем. Для решения ИУ (27) при условии 4(/.2) = 0, (1.2) е © используем ме- 
тод нелинейных граничных ИУ типа Гаммерштейна [1, 2], позволяющий одновременно опреде- 
лить область контакта, контактное давление и нормальное перемещение упругого материала вне 
области контакта. 


Введём обозначения М =(г,2), М=(х,у) и предположим, что область контакта цели- 
ком содержится в прямоугольнике $ = {г -а|<Ь, [2|5 с} ‚ а>Ь>с, который не выходит на 


ребро клина. 

Уравнение (27) дополним условием неотрицательности контактного давления в области 
контакта, а также условиями отсутствия контакта и обращения в нуль давления в дополнительной 
области $ \ © ‚ записав их все в виде системы: 


[к(м,м)а(м)ам =а(м), а(м)>0, Ме®, 


(30) 
[к(м,м)а(м)ам >а(м), а(М)=0, Меб\ о. 
5 
Введём нелинейные операторы 
р’ (М) =зир{р(М),0}, р’(М)=шмЕ{р(М),0}. (31) 
Идея метода заключается в представлении искомого давления в форме 
9=9(М)=9’(М)+а (М) (32) 


с целью автоматического удовлетворения интегрального неравенства (33) в ходе решения нели- 
нейного операторного уравнения типа Гаммерштейна 


9р=0 (МеО), @р=р’+Кр' -4, (33) 


где р= (М), р’ = Р(М), а= &М), 
Кр” = [К (м,М) р" (М)ам. (34) 


Эквивалентность системы (30) и уравнения (33) устанавливается в теореме, аналогичной извест- 
ным теоремам [1, 2]. При численном решении уравнения (33) применим модифицированный ме- 
тод Ньютона [1]. 

Введём следующие безразмерные обозначения: 


г-а=г'6, х-а=х'6, 2=2'6, у=у'Ь, 6б=0'Ь, (35) 
А=Ь/(2к,), в=Ь/(2в,), Л=а/Ь, к=ВЬ, Е=СИЬ, (36) 
9'("',7')=9(г,2)/(2п6), Р’=Р/(2п0Ь?), 5'>5, 99. 87 
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Штрихи далее опускаем. Параметр Л характеризует близость областей контакта к ребру клина, 
параметр к характеризует удалённость штампов друг от друга. В расчётах полагали 
б=А=В=Е= 1, у= 0,5. В таблице даны значения вдавливающей силы Р для разных углов клина 
2а при изменяющихся значениях Л, ки р. При движении штампа к ребру клина (р<0) сила мень- 
ше, чем для случая удаления от ребра (>0). С уменьшением значения ^ или угла клина величина 
Р падает. Эти факты объясняются приближением областей контакта к грани, где отсутствуют 
напряжения. Известно [3], что вдали от ребра клина величина ДА) практически не зависит от ко- 
эффициента трения |. Вблизи ребра клина это уже не так. Чем больше осадка, тем больше сила. 
Большая вдавливающая сила требуется при удалении штампов друг от друга (параметр к возрас- 
тает). Значение к = со соответствует случаю одного штампа. Взаимовлияние штампов существен- 
но проявляется при уменьшении параметра к (36). 
Значения вдавливающей силы Р 






































к 1,2 2 3 со 1,2 2 3 со 

|1] 2а = 60°; Л = 1,2 2а = 60°; Л=2 

0,2 0,071 0,079 0,083 0,102 0,092 0,098 0,105 0,134 

0 0,064 0,070 0,075 0,092 0,082 0,089 0,096 0,123 

—0,2 0,058 0,065 0,069 0,084 0,076 0,082 0,088 0,113 
2а = 90°; Л = 1,2 2а = 90°; лЛ=2 

0,2 0,172 0,185 0,195 0,221 0,193 0,206 0,215 0,244 

0 0,157 0,170 0,180 0,205 0,179 0,192 0,202 0,233 

—0,2 0,145 0,157 0,166 0,191 0,168 0,180 0,190 0,222 





























Заключение. Решена новая пространственная контактная задача с неизвестной областью кон- 
такта для упругого клина о взаимодействии двух одинаковых эллиптических в плане штампов с 
учётом трения при одной свободной от напряжения грани клина. При использовании метода не- 
линейных граничных ИУ с учётом симметрии задач определены области контакта, давления в 
этих областях, вдавливающие силы. Сделаны расчёты при разных значениях угла раствора клина 
и коэффициента трения. 
Работа поддержана грантами РФФИ 12-01-00065, 12-01-00991. 
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ЗО СОМТАСТ РКОВЕЕМ ОМ ТМТЕВАСТТОМ ВЕТ\ММЕЕМ РЕЕХТВЕЕ ММЕОСЕ АМО ТМ/О 
РУМСНЕ$ УМОЕК ЕВТСТТОМ 


О. А. Ро2Ваг$Юу, А. $. Ро[уаКом 
(Роп За Тесйптка! Упмег®Ку) 


Тре 5райа!/ сотасЕ ргоМет юг {те ПехЫе иеаде /5 туезНдаеа т {те диаз!-$айс юЮгтШайоп ипеп мо 1аепёса! 
еШрёс пога рипсвез аге збуттейтсайу /тргеззе4 тю опе Гасе оЁ {Те иеаде ипаЕег тсйоп. Тре о ег иеацде Гасе 
[5 $е55-ее. Тре сотасЁ дотат 15 ассереа ипкпоит. Тпе тойоп о! {те рипсвез регрепасишаг №0 1е едде о! 
{е иеадде 15 сопу/дегеа, [п #5 сазе те рипсйез сап тоуе ейпПЕГ {оюиагаб {те еаде ог Гот те еабе. Изтд те 
Кпоитп Сгееп Гипсвоп Юг 1е гее-атепзюпа! иеаде, {Те сощасЁ ргоМет /5 гедисей ю {те тедга/ едиайоп 
иИёИ гезресЁ ю те сопёасЕ ргеззиге. баапоу5 те{то4 1/5 изеа Юг оМпда ве теодга! едиайоп иисй регтй5 ю 
Ч&егтте {Те сотасЕ аотат апа {те сожасЕ ргеззиге хти!апеоиз. ТВе питегса! апа/у/5 [5 таде юг аИТегепЕ 
уа/иез о! ве своп соейсепЕ иеаде апдЕ, апа те @апсе Белеет те рипсвез. 

Кеуигогаб&: еогу оРе$Нсйу, сопёасЕ ргоМет, игеаде, своп, ба/апоу5 тетод. 
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